Dynamics of polymeric manifolds in melts: the Hartree approximation 
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The Martin-Siggia-Rose functional technique and the selfconsistent Hartree approximation is 
applied to the dynamics of a D-dimensional manifold in a melt of similar manifolds. The generalized 
Rouse equation is derived and its static and dynamic properties are studied. The static upper 
critical dimension, d uc = 2D/(2 — D), discriminates between Gaussian (or screened) and non- 
Gaussian regimes, whereas its dynamical counterpart, d uc = 2d uc , discriminates between Rouse- 
and renormalized-Rouse behavior. The Rouse modes correlation function in a stretched exponential 
form and the dynamical exponents are calculated explicitly. The special case of linear chains D = 1 
shows agreement with Monte-Carlo-Simulations. 
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Much attention has been paid recently to the theory of 
the dynamical behavior of polymers (or generally speak- 
ing polymeric manifolds) and flux-lines in a quenched 
disordered random medium ]3] or in a melt . Tech- 
nically either the projection formalism and the mode- 
coupling approximation J0|| or the selfconsistent Hartree 
approximation (HA) have been used for the deriva- 
tion of the equations of motion. 

In this paper, the HA is used to investigate the static 
and mainly the dynamic properties of a polymeric D - 
dimensional manifold (or a fractal) in the melt of the 
same manifolds. The generalized Rouse equation (GRE), 
which we derived, reproduce in the static limit the screen- 
ing and saturation of D - dimensional manifolds || in 
a different way. The whole dynamical consideration re- 
sults in a subdiffusive behavior and exponents, which are 
confirmed for the 3-dimensional melt by MC and MD 
simulations |^-^|. We should stress, that the manifolds 
in our consideration are crossable, so that entanglements 
cannot occur and reptation dynamics is not considered. 
The dynamics considered here corresponds to chains be- 
low the critical molecular weight: entanglements are not 
of importance for the dynamic behavior [0] . We describe 
the manifolds below only in terms of connectivity and 
excluded volume. The connectivity defines the D - di- 
mensional subspace which is embedded in the Euclidean 
space of d dimension. The model we have chosen allows 
the interpolation between linear polymer chains, which 
correspond to D = 1, and tethered membranes (D = 2). 
By analytic continuation to rational numbers of the spec- 
tral dimension statements on polymeric fractals can be 
made. In a series of papers |3|, we have considered the 
different regimes in static scaling. In the present letter 
we place emphasis on the dynamics. We will show be- 
low, that a new dynamical regime for the motion of the 
manifold segments appears. 

Let us start with a melt of D-dimensional manifolds 
in a d-dimensional space. The test manifold is repre- 
sented by the d-dimensional vector R(x, t) with the D- 
dimensional vector x of the internal coordinates. In the 
same way the manifolds of the surrounding matrix are 



specified by r^- p \x, t) (p = 1,2, . . . , M). We have chosen 
the notation in such a way, that the boldfaced characters 
describe the external degrees of freedom in Euclidian d 
- dimensional space, whereas the arrow hatted vectors 
correspond to the internal D - dimensional space. The 
model of the melt of M (monodisperse) tethered man- 
ifolds used in the following is based on the generalized 
Edwards Hamiltonian, 
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In this melt an additional (test) manifold is immersed 
which is described by the variables R(x) . The number of 
monomers along one side of the manifold is N and limits 
the x - integration. 

The corresponding Langevin equations in Cartesian 
components j for the test chain has the form 
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and similarly for all other polymers in the melt 
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where £o is the bare friction coefficient, e = Td/l 2 the 
elastic modulus with the Kuhn segment length /, V{- ■ ■ ) 
the excluded volume interaction function, denotes a 
-D-dimensional Laplacian and the random forces fj and 
fj have the standard Gaussian distribution. 

We find it more convenient to reformulate the Langevin 
problem (2,3) in the MSR-functional integral representa- 
tion [|. This representation is especially useful for per- 
forming transformations to collective variables or integra- 
tion over a subset of variables. In our case we introduce 
the matrix density p(r,t) and the response field density 
7r(r,i) 

M . 

p(r,t) =Y, / d D xS(r-r ( - p \cc,t)) (4) 

M d 
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In ref. |[o| the first systematic expansion of the effective 
action in the MSR-functional integral in terms of p and 
7r was given. 

The aim now is to integrate over the matrix variables 
(4,5). To do this, we make the expansion of the effec- 
tive Action up to the 2-nd order with respect to p and 
7T, which corresponds to the random phase approxima- 
tion (RPA). After performing the (Gaussian) functional 
integration all information about the matrix is comprised 
into the RPA correlation Soo (k, t) and response Soi (k, t) 
functions p0| . 

The resulting Action includes still the test manifold 
variables in a highly non-linear way. In order to han- 
dle it we use the Hartree-type approximation and also 
take into account the fluctuation-dissipation theorem for 
both, the test manifold and the matrix variables. This 
strategy leads (see for details in Q) to the following GRE 

d 

to-T^Rjix, t) 

+ j d D x' dt'T(cc,x';t-t')-^pR j (x',t') 

d D x'fl(x, g)Rj{x', t) = Fj{x, t) (6) 
with the memory function 

x^(k;f,f';t)5 00 (k,t) (7) 

(where F(k; x, x';t) is the test manifold density-density 
correlator), the effective static elastic susceptibility 
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x F st (k; 2 , x 1 ) - S(x - x 1 ) I d D x"F st (k; £, x")) (8) 



and the random force has the correlator 

{T^x, t)T 3 {x' , t')) = 2T8 ij k S(x - x')6(t - t') 

+ 6{t-t')T(x,x";t-t') (9) 

In eq. (8) the effective potential 

V(k) = V(k) [1 - V(k)S st (k)/T] (10) 



gains the standard screened form 



V(k) = V(k) l + V(k)F^(k)/T 



(11) 



(where F st (k) is the free system correlator) if the stan- 
dard RPA-result is used for the static correlator S st (k). 

It is an important point that we treat on an equal 
footing both, the static and dynamic parts of the GRE 
(6). Let us start from the static behavior of eqs. (6)-([TTl). 
The static limit of these equations for the Rouse mode 
correlator, C(p) = (R(p) • TL(—p)), yields the Dyson-like 
form 



C(P) 
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where M = N D is the total number of monomers in the 
manifold and the "self-energy" is given by 
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The test manifold static correlator in eq. (|12|) is param- 
eterized by the wandering exponent £ in such a way 



{2ix) d l + V(k)F^\k)/T 
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In its turn the exponent C is determined by the correlator 
C(p): 
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The system of eqs. (|L 2|) - (|l 5|) should be analyzed self- 
consistently. Straightforward calculations yield 



C,(d+d° { +2)-D 



(16) 



with the Gaussian fractal dimension d® = 2D/(2 — D) 
[p|,^2[ and c%, C2 are some constants depending on d and 
D. Physically, the condition for the exponent ( comes 
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from the balance between the entropic and the interac- 
tion terms in the denominator of eq. (|l2]). The anal- 
ysis shows that the only way to satisfy the eqs. (|l2|)- 
( |i~5| ) is to impose the exponents in eq. (|l6|) the condi- 
tion: C(d + 2) = ((d + d° + 2) - D > 2. This holds 
at d > d uc = 2D/(2 — D) and the only solution is 
the Gaussian one with £ = £o = (2 — D)/2. Besides 
that, the necessary condition d® < d immediately implies 
D < D s = 2d/(2 + d). The upper critical dimension d uc 
in a melt and the spectral critical dimension D s was dis- 
cussed first in j^]. At d < d uc the interaction term in 
eq. ( |l2| ) overwhelms, |S(p)| ^ (2irp/N) 2 , and the system 
become unstable. The manifold is saturated in a melt, 
i.e. it loses its fractal nature and becomes compact ||. 

We now consider the dynamics at d > d uc - There are 
two dynamic exponents, z and w. The exponent z mea- 
sures the time dependence of a monomer displacement, 



Q(t)= / d D p / ^L[l- e st ]C(f,s)oct 2z 
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and the exponent w measures the same for the center 
of mass 



Qc.m.(i) = lim / 
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In eqs. (0)-© C(p,s) = (|R(p,s)| 2 ) is the Rouse- 
Laplace component of the correlator C(x, t). The formal 
solution of eqs. (6)- (8) for C(p, s) is given by 



C(p, s) = 



Cst(p) 



6 £ +Afr(p,s) 



(19) 



where C st (p) = l 2 {N/2-Kp) 2 jM and r(p, s) is the Rouse- 
Laplace transformation of the memory function (7). 

In the RPA the matrix density correlator Sqq (k, t) is 
well approximated by 

Soo(M = S st (k) 

exp{-fc 2 £ coop (k)i} , {kl) d °N«.l 

x< «*{-«(*)*"} ■ («H">i (20) 

where the cooperative diffusion coefficient -D coop (k) = 
pV(k)/£ , r = Co' 2 /rd and z = (2- D)/4 is the Gaus- 
sian z-exponent. The corresponding Ansatz for the test 
manifold yields 



F(k;x;t) =F at (k,x) 
e?q>{-k 2 D G t} 
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(i) By making use of eqs. (20)-(21) and eq. (7) in the 
limit {kl) d °fN > 1 one can derive the result 



A/T(p, s) — const. 
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with 



/3 = z Q (d - d uc + 2) 



(22) 



(23) 



In the derivation of eq. ( P2| ) we also used the assumption 

l 2 <C Q{t) < l 2 J\f 2/d ° . The condition [3 < 1 immediately 
defines the dynamical upper critical dimension 
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i.e. the dimension above which the manifold has the 
simple Rouse behavior, at d = d uc one can call it the 
marginal Rouse behavior and only at d < d uc the dy- 
namic exponents z and w are renormalized. 

At d < d uc eqs. (|l^), after inverse Laplace transforma- 
tion, yields 



oo 



T(m/3+l) (25) 



where 4 = [(3\V(k = O)! 2 ^^ = 0)/l d+2 ] and T(x) is 
the gamma-function. The eq. (25) is very close to the 
stretched exponential form found by the MC- and MD- 
simulations ||. The eq. (25) was actually calculated in 
the limit p — > 0, and we can use it first of all to compar- 
ison with simulation results on the center of mass mean 
square displacement. By using eq. (25) in eq. (18) we 
obtain 
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where V = l 2 eA/T(/3 + 1) and 

w — — z (d - d uc + 2) 



(26) 



(27) 



By the same way from eq. ( |17j ) we can derive: 

z = Zq/3 = z 2 {d - d uc + 2) (28) 

In ref. [3] this renormalized dynamics was formally used 
as a projected dynamics for the polymer mode coupling 
approximation (PMCA) and eventually leads to the GRE 
which, as the authors claim, can describe the entangled 
dynamics. We will argue in an extended paper [p"l| that 
this is a result of misinterpretation of the GRE. 

(ii) If we assume that the main contribution to the inte- 
gral (7) comes from the small wave vectors, {kl) d fN <C 1, 
then we arrive at 



where Da is the self-diffusion coefficient. 



A/T(p, s) oc s ( rf -<*-)/2 



(29) 
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Since d > d uc , the simple Rouse behavior in the small 
wave vector regime does not change. 

At the large displacement regime, R G <§C Q c .m{t), one 
should expect a simple diffusive behavior: 

Qc.m.(t) = dD G t (30) 

with 

D G = T/Mfa + NT(p = 0; s = 0)] (31) 

Now the problem is how Dq depends from Af? One 
can assume that for this case only the small wave vec- 
tors, (kl) d <AT < 1, are relevant, i.e. the dynamics of 
the matrix is driven by the cooperative diffusion coeffi- 
cient -Dcoop and the dynamics of the test manifold by the 
self- diffusion one Dq (see eq. (20-21)). Since in any way 
-Dcoop ^> Dq, the calculation yields 

ATr(p = 0; S = 0) oc [^coopl-V 1 "^' (32) 




T-R log(T) 



FIG. 1. A schematic plot of Q C m(t) for the simple Rouse 
(dashed line) and the renormalized Rouse (solid line) dynam- 
ics 



length N — 200 at the relatively short time Q C m(t) oc t w 
with w = 0.8 (instead of w = 1) in || and w = 0.71 
in 0| . This deviation from the simple Rouse regime also 
occurs for short chains (N < N e ) which clearly are not 
entangled ||. 

The best test of the renormalized Rouse dynamics pre- 
dictions would be the simulation of rather long crossable 
(to avoid reptation) chains but still with an excluded vol- 
ume interaction. In a recent MC-simulation Jl3|] the stat- 
ics and dynamics of such melts have been studied. Un- 
fortunately in |l3] | the plot Q C m(t) is not given explicitly, 
i.e. it stays unclear from this simulation how the mode 
p — > is renormalized. 

In summary, we have shown that by using MSR- 
functional technique and the Hartree approximation the 
GRE for a D-dimensional manifold in the melt of similar 
manifolds can be derived. In this equation the static and 
dynamic parts are treated on an equal footing. Besides 
the static upper critical dimension, d uc = 2D/(2 — D), 
its dynamical counterpart, d uc = 2d uc was found (such 
that at d uc < d < d uc the manifold is Gaussian but 
renormalized-Rouseian) . We have calculated the dynam- 
ical exponents, w and z, and have explained some novel 
computer simulation findings. 

We have benefited from discussions with J. Baschnagel, 
K. Binder, K. Kremer, W. Paul and T.B. Liverpool. We 
are greatful to the Sonderforschungsbereich 262 and the 
Bundesminister fur Bildung und Forschung for finacial 
support. 



But .Dcoop = O(N ) and <f f < d, then T(p = 0; s = 
0) -> at M -> oo. As a result Dq = T/N£ , i.e. the 
simple Rouse result does not change. 

In Fig. 1 we have summarized the overall schematic 
behavior for Q cm (t). At the relatively short times, tq < 
t -C is, and displacements, I 2 < Q cra {t) <C Rq, the test 
chain dynamics is mainly determined by the fluctuations 
from the interval (kl) d fj\f ^> 1. As a result the renor- 
malized Rouseian behavior dominates and e.g. for the 
melt of polymer chains [D = 1) in the 3 - dim. space 
the exponent w = 3/4 = 0.75. In the opposite limit, 
Tfl <C t and R G ^ Q C m{t), the long wavelength fluctua- 
tions {klffAf < 1 are relevant and the melt almost does 
not influence the test chain: The simple Rouse regime is 
recovered. 

MC-simulations of the bond-fluctuation model ^| as 
well as the MD-simulations of the athermal melt have 
been undertaken. Recently also the static and dynamic 
properties of a realistic polyethylene melt have been stud- 
ied Q . Both in MC and MD simulations a slowed down 
motion at intermediate times for the center of mass is 
clearly observable. It was found e.g. that for the chain 
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